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1. Introduction 

In this paper we study solutions of the following model equation for the time-dependent 
vector field u(x, t) — (u\{x, t), . . . , u n (x, t)) on E™ x (0, T) 

du 1 1 

— +au\7u + -(1 - a)V|u| 2 + -(divu)u = Au + xVdivu, (1.1) 

where a € (0,1) and n > are given parameters. The equation l|l.lfl is of interest for 
various reasons. For example, it has the same scaling properties and the same energy 
estimate as the Navier-Stokes equation (NSE): If u(x, t) is a solution of p. 1(1 then also 
Au(Ax, \ 2 t) is a solution for A > and, for sufficiently regular solutions with a suitable 
decay at infinity, we have 

\\u{ Xl t)\ 2 dx = ( \u{x,t')\ 2 dx (1.2) 



/ ( [|Vu(x, s)\ 2 + x{d\Yu{x,s)) 2 ] dxds. 



Heuristically, solutions of (|1.1|) should converge to the solutions of the NSE as h — ► oo. 
Similar penalization schemes have been used in numerical algorithms for solution of 
NSE, see, e.g., 0. 

In dimension n = 2 equation i|l.lf> is "critical" (i.e. the controlled quantities 
are invariant under the scaling symmetries of the equation) and hence it is natural to 
expect that the full regularity of solutions with finite energy can be proved by standard 
methods. 

In this paper we shall concentrate on the super-critical case n > 3. It is natural 
to expect that the theory of Leray's weak solutions applies in this case. Moreover, it 
is likely that for n = 3 the partial regularity results in the spirit of Scheffer an d 
Caffarelli-Kohn-Nirenberg |3] can be proved here as well. We note that for a = 1/2 
the non-linear part in (|l.lf> can be written in the divergence form and consequently 
one can directly apply the known regularity theory for the NSE in that case. However, 
most questions regarding full regularity of solutions to Hl.lfl in the case n > 3 appear 
to be open. 

Our aim here is to investigate the problem of finite-time blow-up for a special 
class of solutions to 1(1. ljl . We study solutions given by 

u(x, t) = — w(r, t)x , (1-3) 

where r = \x\, and v(r,t) is a scalar function. Such vector fields, usually called radial 
vector fields, are not divergence free unless v = and hence the relevance of such 
solutions for the theory of the NSE may be limited. Nevertheless, the behaviour 
of these solutions provides an interesting insight into various scenarios of singularity 
formation. 

Using the radial vector field ansatz and substituting in the equation 1(1. ljl we 
obtain 

/ n _i_ \ \ 

Vt = (1 + x) v rr H v r + 3rvv r + (n + 2)v 2 , (1.4) 



where subscripts denote corresponding partial derivatives. Replacing v(r,t) by 
(1 + x)v{r, (1 + x)t) we see that, when studying the radial solutions, one can assume 
k = without loss of generality. 
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Our first result is that in dimension n < 4 the solutions to (|1.4|l do not exhibit 
blow-up if there exists C > such that the initial condition vOr, 0) = Uo(r) satisfies 
- C < v (r) < C*(l + r)-(™+ 2 )/ 3 when n < 4, and (1.5) 

-(- 

r 2 I 3 



-C < v (r) < — (- log r + C) when n = 4 . (1.6) 



On the other hand, when v(r,0) — vo(r) — c > (where c is a constant), then 
v(r, t) — v(t) solves 

^L = (n + 2)v 2 , v(0) = c, 

and the solution blows up at time t = l/(c(n + 2)). Therefore some control of Vo(r) 
at infinity is necessary to prevent formation of singularities. 

The proof, that conditions (|1.5|) - (|1.6[) are sufficient for preventing blow-up, is 
based on an analysis of steady-state solutions to the equation (|1.4|) . The steady states 
can be analyzed more or less completely since the equation 

71+1 

v" + v' + 3rvv' + On + 2)v 2 = , (1.7) 

r 

can be transformed to an autonomous two-dimensional dynamical system. We briefly 
outline the behaviour of the steady-state solutions: Equation (|1.7fl has a solution 
V: [0, oo) — * (0, oo ) with V(0) = 1, V^'(O) = and the following asymptotics at infinity 

V(r) - r~ (n+2)/3 , when 1 < n < 4 
VOr) ~ r~ 2 ^ logr + C^j , when n = 4 

VOr) ~ r~ 2 , when n > 4 . 
Using the scaling symmetry we obtain a one-parameter family of solutions 
v x (r) = X 2 V(Xr). 

There are also other interesting steady-state solutions. It turns out that even for 
radial solutions, weak solutions of l|l.l|) can exhibit the following non-trivial behaviour: 

(i) formation of singularities with a different rate of blow-up than suggested by scaling, 

(ii) violation of local energy inequality, 

(iii) significant non- uniqueness. 

We present more specific discussion of these phenomena in Section [21 

The second group of the results, we shall discuss, concerns the blow-up behaviour 
of solutions to the equation (|1.4(l in dimensions n > 4. In this case our results are 
based on combination of analytical arguments and numerical calculations. We will 
present strong evidence that for n > 4 and suitable compactly supported initial data 
there exist solutions of (|1.4|) that form a singularity in finite time. 

A natural class of singularities for the equation (|1.4|) are self-similar singularities 
of the type 

1 ( r \ 

v(r,t) 



2n(T-t) \y/2 K (T-t) J 

where k > 0, T > are parameters and w is a function defined on [0, oo). The 
equation for w is 

71 + 1 

w" H w' — nrw 1 + 3rww' + (n + 2)w 2 — kw = , (1-8) 
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together with the natural boundary conditions 

w(0)=a>0, w'(0) = 0, (1.9) 

w(r) ~ r~ 2 , as r — > oo. (1-10) 

The scaling symmetry (w, k) — > (X 2 w(Xr), X 2 k) allows us to fix k and vary only a. 
There are no non-trivial solutions of (|1.8fl - l|1.10|) f° r 71 — 4- We conjecture that for 
n > 4 and a fixed k, the number of non-trivial solutions is determined by the number 
of roots of the solution to a linearization of (| 1 - H|l around its trivial equilibrium 

2k 

w = . 

n + 2 

The relevant linear problem is 

„ (n + 1 n - 4\ , , 

z"+[ nr -)z' + 2nz = 1.11) 

\ r n + 2 J y ' 

z(0) = 1 , z'(0) = 0. 

The solution of this problem can be written explicitly in terms of a confluent hyper- 
geometric function: 

. . „ r ( n + 2 n + 2 r 2 

z(r) = M , , k— 

y ' \ n-4' 2 ' 2 

We refer the reader to for basic properties of hypergeometric functions. It turns 
out that z(r) is a polynomial whenever (n + 2)/(n — 4) is an integer. The number 
of zeros m = m(n) of z(r) in (0, oo) is the smallest integer greater or equal to 
(n + 2)/(n — 4). We recall that we assume n > 4 at this point. We conjecture 
that, for a fixed K and n > 4, the number of solutions of (|1.8|) - (|1.1U|) is m(n) — 2. 

There are certain similarities between the behaviour of solutions of (|1.4|) and 
solutions of a system arising in mathematical biology studied in [P . 

In the last section we study some related problems concerning a widely studied 
semi-linear heat equation. There appear to be striking similarities between the 
formation of singularities in the equation 11.4JI for n > 4 and similar behaviour of 
non-negative radial solutions of 

— = Av + v 2a+1 , inR n x(0,t) (1.12) 



for n > 10 and a > a c (n) = 2/(n — 4 — 2y/n — 1). For the significance of the critical 
exponent a c (n) see, for example, [HlC]- 

Self-similar singular solutions of l|1.12|l are given by 



v(x,t) = (2k(T -t))- 1/a 



w 



y/2 K (T-t) 



where k,T > are parameters and w is a function on [0, oo). The function w is a 
solution of the boundary value problem 

" " -\ w > - K (™ +rw >\ +w 2-+i =0 (1.13) 



r v (j 



w(0) — a , w'(0) = , and w(r) ~ r ' CT as r — > oo (1-14) 
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The existence of non-trivial solutions to (|1.13|1 - (|1.14|1 depends on n and a in the 
following way. If 

2 _ 

a < , there are no non-trivial solutions, see HI, 

n — 2 

2 2 . _ 

< a < ^=^= , there are infinitely many solutions, see 1121 . 

n — 2 n — 4 — 2^ - 1 

2 3 _ 

< (j < , there exists at least one non-trivial solution, see 19]. 

ra- 2 n - 10 > uj 

One of the open problems for 14|t is to determine the exact range of 

parameters for which the boundary-value problem has a non-trivial solution. In 
Section [3] we present strong evidence that the sufficient condition of Lcpin (UJ) 
a < 3/(7i — 10) is also necessary for the existence of non-trivial positive solutions 
of (031 - 1011) . 



2. Solutions of the model equation 

2.1. Phase portrait 

In this section we analyze the steady-state solutions of the equation 

v" + — -v' + 3rW + (n + 2)v 2 = . (2.1) 
r 

The invariance of solutions to l|2.1|l under the scaling v(r) — > A 2 w(Ar) suggests the 
change of variables: 

v = r~ 2 w , r = e s , 

which transforms l|2.1|l into an autonomous equation 

w" + 3W - (4 - n)w 2 + 2(ti - 2)w - (4 - n)w' = . (2.2) 

With a slight abuse of notation the prime ' now denotes differentiation with respect 
to the new independent variable s. 

We are interested in the phase portrait of the vector field in R 2 defined by (|2.2|) . 
Global properties are best studied in a suitable compactification of M 2 . It turns out 
that the transformation of variables 

/ \ ,,, „ , , , t&nib(s) 

w(s )=tan(0(» , w'(s) = ^ \ , 

cos^ <p[s) 

leads to a compactification which works well in the case at hand. In the new variables 
the equation 1)2.2(1 becomes 

' -PiM), (2.3) 



ds cos <f> cos ip 
dtp 1 



Q{<t>,1>), 



ds cos 4> cos tp 
where 

P(4>,tp) = cos 4> sin ijj , 

Q(4>, ">P) = — sin ip cos ip sin <^>(2 sin?/; -I- 3 cos?/>) + 2(n — 2) sin ^ cos 2 (/icos 3 i 
+ (4 — n) cos 2 ip cos </>(sin 2 ^ cos ip + sin -0) . 
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Table 1. Equilibria of the system 12.41 in the region (~7r/2, 7r/2] x (—it/2, it/2] 

Equilibrium (<j>, ip) Linearization 

eigenvalues, eigenvectors 

ev- (01,^0 = (0,0) Al=2 ' ( 2 ) ' A 2 = -(™- 2 )' ( n\ 

e 2 : (02, ^a) = (arctan 2( ^T 4 2) , 0) Ai,2 > if n < 4 , Ai,2 < if n > 4, 

e 3 : (03,^3) = (f,0) Ai = -3, ( J ) , A 2 = 0, f ^ 

e 4 : (^4,^4) = arctan (-§)) Ai = 1 , A 2 = ^ cos 3 ?/>4 > 1 

e 5 : (05,^5) = (f,f) Ai = -1, ( J J , A 2 = 2, f M 



Hence the integral curves of the vector field 1)2. 3J1 on the torus T 2 = M 2 /27rZ 2 ar e 
defined by 

^=P(0,VO, ^ = Q(«^). (2-4) 
as as 

Due to the periodicity we have 

P(cj) + 7T, tf) = P(<P, iP + tt) = -P(cj>, i>) 

Q(<t> + 7T, V) = Q(0, V + t) = ■ 

Therefore it is sufficient to analyze the flow in the region \<f>\ < tt/2, \ip\ < tt/2. 
We proceed with a description of equilibria and important heteroclinic orbits in the 
region (— tt/2, tt/2] x (— w/2, tt/2]. Both, the orbits and the equilibria are computed 
in a fully rigorous way by standard methods. All the other equilibria are obtained by 
shifts along the coordinate axes by kit. From the information about the equilibria and 
the heteroclinic orbits we determine the full phase portrait of the system Q2.4[l , The 
phase portraits differ for different values of the dimension and we sketch the three 
different cases n < 4, n — 4 and n > 4. Important heteroclinic connections in a 
typical phase portrait when n > 4 are depicted in Figure ^a), where n = 5 was used 
in numerical computations. Similarly Figure^b) depicts the phase portrait for n = 4 
and Figure ^c) for n < 4 (computed for n = 3). The curves denoted by a, b, c, d, / 
are important in our considerations and their meaning is explained below. 

2.2. Specific solutions 

Analysis in the previous section enables us to show existence of some important 
solutions of 12. f fl . 

Orbit (a) The heteroclinic orbit (a) in Figure Ufa) gives immediately the solution 

v : [0,oo) -> R, v > 0, 

v(0) = 1, 1/(0) = 0, 
v(r) ~ r~ 2 , as r — > oo . 

Orbit (b) The orbit (b) in the phase portrait for n > 4 (Figure ^a)) corresponds to 
the solution v : [0, oo) — > K, v > 0, 

v(r) ~ r ~ , as r — > 0+ , 
v(r) ~ r~ 2 , as r — ► oo . 




The behaviour as r — ► 0+ is determined from the slope of the center manifold at 
the equilibrium point e%. If (<fi,ip) approaches along (b) we have the following 
asymptotic expansion 

^=^(0-|)+a 2 (^-|) 2 + ... + a fc (0-|) fc ... (2.5) 

The series is not convergent as is typical for similar situations involving center 
manifolds. The coefficients cik can be calculated but we will not need their exact 
values for the subsequent analysis. The expansion (|2.5(l gives w(s) — > oo as s — > — oo 
and 

w'(s) = -'^-w(s) +b + b-1.1v- 1 (a) + ... + b- k w- k (s) ... , (2.6) 

which again is an asymptotic expansion, not necessarily convergent. From Ij2.6|l we 
can see the described behaviour of v(r). 

Orbit (c) A similar analysis of the orbit (c) in Figure da) gives the solution 

v : [0,oo) -► R 

v(r) ~ r ~ , as r — > 0+ , 
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v(r) ~ r~ 2 , as r — * oo . 

Orbit (d) Exploiting again the known slope of the center manifold at e-$ one easily 
sees that the heteroclinic orbit (d) in the phase portrait for n < 4 corresponds to the 
solution v : [0, oo) — > R 

u(0) = 1, v'(0) = 0, 

v(r) ~ r 3 , as r — ► oo . 

Orbit (f ) In the case n = 4 we have an explicit equation for the solution corresponding 
to the orbit (f) 

1,4, 1 w 2 

——w + — log ■ 



3 9 b l-3u//4 2 ' 
which leads to the solution v : [0, oo) — ► R 
u(0) = l, i/(0) = 0, 

u(r) ~ r~ 2 I — log r + C ) , as r — * oo . 



The existence of solutions (d) and (f) implies the following result: 

Theorem 2.1 Assume that uq(x) = — vq(\x\)x/\x\ is a smooth radial vector field on 
R™ and assume that one of the following conditions is satisfied 

(i) n = 2, 3 and 

-C < «o(M) <C(l + |a;|r ( ™- 1)/3 , for some C > 0, 

(ii) n = 4 and 



(m) n > 5 and 



^ <r n 4/31og(l + |xl)+C 

-C < u (f )< > for some C > 0, 

1 + x 



-C < fo(|x|) < 7^ ttt-T) for some j < 1 and C > 0. 

(n-4) |x| 



TTien </ie equation with x = /ias a global bounded solution with the initial 

condition u . 

Proof: The proof of the above theorem follows from the analysis of the two- 
dimensional system discussed above. We use the solutions constructed in this analysis 
as barriers in the equation (|1.4() . 

Remark 2.1 As we have noted in the introduction, when studying radial solutions to 
the equation 07j e can consider x = without loss of generality. 

2. 3. Singular steady states for n > 5 

We consider the following boundary- value problem for radial vector fields u : R" — » R": 

- Au + -uVu + iv— + -udivu = , in B± (2.7) 
2 2 2 2 v ' 

u(x) = -bx, on &Bi , (2.8) 
where B x = {x £ R™ | |x| < 1} denotes the unit ball and 6 £ R. 
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Definition 2.1 We say that u is a weak solution of [2.7jl if it belongs to W ' 2 (Bi), 
satisfies the equation in the sense of distributions, and the boundary condition is 
satisfied in the sense of traces. 

Furthermore, we say that u is a suitable weak solution of \2. Tjj if it is a weak solution 
and satisfies a local version of the energy inequality 



Vu| 2 0-iM 2 A0-i( M .V0)|u| 2 

I Bi L 1 L 

for each smooth <j) > compactly supported in B\ 



<0, (2.9) 



It is easy to see that radial weak solutions are smooth away from the origin. For 
solving l(2~7jl - lf2~5|l we can use orbits (a) and (c) discussed in the previous section. 
The orbit (a) gives a smooth solution when < b < 2(n — 2)/(n — 4). When 
b > 2(n — 2)/(n — 4), problem l|2.7|l - H2.8(l does not have a smooth radial solution. 
In that case a suitable weak solution as defined above can be obtained from the orbit 
(c). It is interesting to note that this solution has a singularity at the origin which 
is not asymptotically self-similar. The self-similar rate of blow-up would be | a; | 1 , 
whereas the actual rate is | a; | ^ x ) / 3 . Moreover, we note that the solution satisfies 
the inequality (|2.9|) . but it does not satisfy the local energy identity, i.e., I|2.9|l with 
equality. 

Apart from the suitable weak solution described above there exist many other 
weak solutions to l|2.7[) when b > 2(n — 2)/(n — 4). These solutions are constructed 
using heteroclinic orbits such as the orbit (c) in the phase portrait for n > 4. However, 
these solutions do not satisfy H2.9|l . 

The above discussion does not cover the case b < 0, which is left to the reader as 
an exercise. 

2.4- Self-similar singular solutions 

In this section we study radial self-similar solutions of (|1.1(1 . i.e. solutions of the form 

where T £ R and n > 0. Assuming (without loss of generality) that x — in the 
equation p.ljl . we obtain the equation 

71+1 

w" H w' — nrw' + irww' + (n + 2)w 2 — nw = . (2-11) 

r 

As is usual in similar situations, it is useful to interpret the equation l|2.11|l as an 
equation of motion for a particle with the unit mass which is moving in a potential 
field, given by the potential V(w) = (n + 2)/3w 3 — kw 2 /2, in the presence of damping 
/j(r, w)w' where fi(r, w) = [n + l)/r + irw — nr. Using this notation we can write 

d 

w" + /i(r, w)w' = — — — V(w) . 

aw 

As in the case of steady-state solutions it is useful to introduce new variables by the 
transformation w(r) — r~ 2 u(r) and r = e s . In the new variables we have the equation 

u" + (n - A)v! + (n - 4)w 2 - 2(n - 2)u + Zuu' = ne 2s u' , (2.12) 

or 

u" + jl(s,u)u' = —^-V{u), (2.13) 
ou 



Singular solutions 10 
where 

p,(s, u) = n - 4 + 3u + «e 2s , V(tt) = m 3 - (n - 2)u 2 . 

Natural boundary conditions for (|2.11|1 in the context of self-similar singularities are 
w(0) = a£l, w'(0) = and (2.14) 
w(r) ~ r~ 2 , as r — > oo . (2.15) 

The condition at infinity together with H2.10|) leads to the blow-up profile 

x 

u(x, T) = —c-, — , for some cEK. 
\x\ 2 

One can easily see that c > in this case. 

The boundary- value problem (|2.11|l . (|2.14(l - (|2.15(l is studied as a non-linear 
eigenvalue problem with parameters a and k considered as unknowns. Due to the 
scaling symmetry (a, k, w) — > (X 2 a, X 2 k, X 2 w(Xr)) we can assume that k = 1 without 
loss of generality. 

Our analysis and numerical computations support the following conjecture 
Conjecture 2.1 Suppose we define 

u(n) = min < k 6 N I k > — — — 
w [ 1 ~ n - 4 

then for a fixed k > anc? n > 4 f/ie boundary-value problem \2.1ty . \2.1$ - \2.1ty has 

u(n) — 2 non-trivial solutions. 

Remark 2.2 (i) Results of Section\£J£ imply that y2.ll)) . \2.14\j - y2.15)) has no non- 
trivial solution for n < 4. 

(ii) Since we impose the boundary conditions \2.11$ and )2.15)) the solutions w = 

and w(r) — 2 [^Z^j i" 2 o,re excluded from the count. 
(hi) The solutions for n — 5 are depicted in Figure^ 

The following simple observation is useful. 

Lemma 2.1 Every solution of the initial value problem )2.11)) and \2. 1$ defined on 
(0, oo) is positive on (0,oo). 

Proof: The lemma follows easily from the interpretation of the problem using the 
equation (|2.13l) . 

We let w a be the solution of the initial-value problem (|2.11|l with the initial 
conditions (|2.14|) and k = 1. Furthermore we define 

R a = sup{i? e (0, oo) | w a > , on (0, R)} . 

By Lemma 12.11 a necessary condition that w a is a solution to the boundary-value 
problem 12.11fl . 12.14fl - 12.15fl is that R a = +oo. Moreover, it is likely that, except 
for perhaps some special values of n, this condition is also sufficient. The reason will 
become apparent later. For n > 4 we let a n = 2/(n+2) and f3 n — 2(n — 2)/(n — 4) and 
we define = j3 n /r 2 . Note that this function satisfies (|2.11|) and also w &n — a n . 

The following formal calculation suggests that w a — > Woo as a — > 0. It is 
convenient to work with the equation 1)2.12(1 . In the coordinates (s,u) the condition 
H2.14fl becomes 

lim er 2s u{s)^a 1 and lim e- 3s (u'{s) - 2u(s)) = . (2.16) 
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Figure 2. Plots of u{r) = —w(r) where w(r) solves 11,81 - 11,101 in dimension 
n = 5 and n = 7/2 (so that the trivial equilibrium is id = 1). The inset shows the 
behaviour of solutions around the equilibrium w = 1. 



-1.5 - 



~-2.5 




We denote by U(a, k, s) the solution of i|2.12[l and l|2.16[l and by (0, s QiK ) the 
maximal interval where s — ► U{a,n, s) is defined. We note that U(e 2s a, k, s) = 
U(a,e~ 2s k,s + s') and will investigate the behaviour of U(a, e~ 2s k,s + s') ass' — > oo. 
The linearization of l|2.12[l at the equilibrium u(s) = [3 n , corresponding to w n in the 
coordinates (r, u) , is defined by the equation 

z" + (n - 4 + ^\ z' + 2(n - 2)z = ne 2s z' . (2.17) 

It is easy to see that the general solution of l|2.17(l is given in terms of confluent 
hypergeometric functions by 

z(s) = de Xls M( ai , h, K^e 2s ) + C 2 e X2S M(a 2 , b 2 , K^e 2s ) , (2.18) 

where a,j = Xj/2, bj = 1 + Xj + n — 4 + 3/3 n and the constants A 2 < — 1 < Ai < are 
roots of the characteristic polynomial for the linear second-order differential operator 
that defines the left hand side of (|2.17|l . The function M is one of the standard 
confluent hyper-geometric functions, see, e.g., For small k > and large s such 
that ne s is controlled, the formula (|2.18() leads to 

U(a, k,s) = j3 n - je Xls M(ai,bi, ^e 2s ) + higher order terms . (2.19) 

We use the approximation H2.19fl to study asymptotic behaviour of U(a, e~ 2s k, s + s') 
for a fixed s and s' — + oo. This is possible since the second term of (|2.19|l becomes 
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_ ry e ^i{s+s )ji(f ( 0l) i^e 2s ) which is small for large s' and fixed s, since Ai < 0. Hence 
we conclude that 

U{ae 2s \ k, s) = U(a, e~ 2s ' k, s + s')=f3 n - e XlS ' je XlS M( ai , b u ^e 2s ) + h.o.t . (2.20) 

This calculation formally shows that w a — > Woo with the rate of convergence 0(a Al / 2 ) 
as a — > oo. This argument can be made rigorous although the details become non- 
trivial. 

For o/(i„, a > 0, we define an index 

i{a) = #{r£ (0, oo) | Woo(r) = a n } . 
The behaviour of i(a) in a neighbourhood of a n is controlled by the linearization of 

da 

1 



the equation l|2.11[) at w &n . Denoting Z — ^r a -\ a =a n we have 



Z" 



+ (3a„ - l)r 



Z' + 2Z = 0, (2.21) 



+(^-x)Z+^Z = 0. (2.23) 



Z(r) = M I — 



r 

with the initial conditions 

Z(0) = 1, Z'(0)=0. (2.22) 

The substitution x = (1 — 3a„)r 2 /2 transforms (|2.21|) into the standard form of the 
confluent hyper-geometric equation 

d 2 Z (n + 2 \ „ n + 2, 

Hence the solution of (|2.21|) with the initial conditions (|2.22f) is found explicitly in the 
form 

n + 2 n + 2 n— 4r 2> 
~4' 2 'n + 2~2~ 

where the function M is one of the fundamental solutions of (|2.23|l (see, for example, 
11 ). Using properties of the function M and recalling that v{n) = min{fc 6 
Z | {n + 2)/(n — 4) < fc} we see that the solution Z has zeros in (0, oo). Before 
stating the following lemma we recall, that R a denotes the first zero of the solution 
w a , and that for a ^ a n , a > i(a) = w a (r) = a n }. 

Lemma 2.2 Assume a\ < ai < cii or < ai < 012 < ot n , R ai , R a2 < 00, and 
i(ai) ^= 1(0.2)- Then there exists a € (0:1,0(2) such that R a — 00 and w a satisfies 

Proof: To prove the lemma we use standard arguments based on the continuity of 
i(a) at points where R a is finite. Some work is required to demonstrate that the 
condition (|2.15() is satisfied, but the arguments are straightforward. 

From the behaviour of the linearized solution at w &n one expects that R a is finite 
near a n , a/a„ and 

lim =( , ^N isodd (2 . 2 4) 

a _> a + [ v{n) + 1 it v(n) is even v ' 

and (2.25) 

lim =( "i n \ + 1 ] { /( n \ isodd (2.26) 
a ^ & ~ I v \ n ) u v \ n ) 1S even 



Singular solutions 



13 



Table 2. Values of a for individual singular solutions tabulated at integral 
dimensions n. The branches of solutions, when continued in n, terminate at 
the dimensions indicated in the last column. There are no solutions with i(a) = 1 
or i(a) > 6. 







n 


= 5 


n 


= 6 


it 


= 7 


n = 8 n = 9 max. n 


<(«) 


= 2 





02631 





09647 





1466 


0.1684 0.17223 n = 10 


i(a) 


= 3 


1 


3830 





2792 





2222 


n= 7 


i(a) 


= 4 





2205 





2505 






n = 6 




= 5 





2940 










n = 11/2 


i(a) 


= 6 





2855 










n = 26/5 



The behaviour of i(a) for large a can be estimated from If2.20fi . The function 
M(ai,b\,x) for ai = Ai/2, b\ = 1 + (Ai +n — 4 + 3/3„)/2 has exactly one root in (0, oo) 
and decays exponentially to negative infinity as i — * oo. This behaviour suggests, 
that the solution w a , for large a, will be "lagging behind" w^,, and will fall back to 
a n for a finite, but large value of r. After that the term rw' will be large and cause 
the solution to reach in finite time. This reasoning leads us to expect 

lim i(a) = 3 . (2.27) 

a — >co 

Similar considerations suggest that 

lim i(a) = 2. (2.28) 

a^0+ 

In this case we use similar heuristic arguments: if a is very close to zero, the solution 
will reach the equilibrium a n at a large "time" r. The term rw'(r) is then already 
significant. After passing through a n , the increase of w will eventually be stopped by 
the non-linear terms and the solution will start returning to a. After it passes through 
a n again, it will reach zero due to the large term rw'(r). Assuming H2.24fl ~ H2.28fl we 
see, that for v(n) odd, i(a) will change from i>(n) to 3 as a moves from a small right 
neighbourhood of a n to infinity and from v(n) + 1 to 2 as a moves from a small left 
neighbourhood of a n to zero. 

Assuming that i{a) changes by 2 we obtain v{n) — 2 solutions. Repeating the 
same argument for v(n) even yields again v{n) — 2 of solutions. 

The above described behaviour and the conjectured number of solutions are fully 
confirmed by numerical computations. In Table |21 we list the solutions for integral 
values of n. 

3. Semi-linear heat equation 

As we remarked in the introduction, there are striking similarities between the singular 
behaviour of solutions to l|1.4fi when n > 4 and the behaviour of singular solutions of 

v t = Av + v 2a+ \ inR n X (h,t 2 ), (3.1) 

for n > 10 and a > a c (n) = 2/(n — 4 — 2\Jn — 1). For example, blow-up solutions of 
13. If) with the blow-up rate (T - i) -1 / " -5 (with some S > for n > 10 and a > a c {n)) 
were constructed in [J]. Such rate of blow-up can be viewed as "slower" than the 
self-similar rate (T — t) _1 / ff . Some authors call this rate "faster" but it seems that 
the term "slower" is more widely used in the present context. It reflects the fact that 
(T — t)~ 1 l' T ~ s becomes infinite "more gradually". We conjecture that the analysis of 
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|Zj can be used for proving existence of slow blow-up solutions for n > 4. In fact, if one 
allows slow decay of initial data at infinity, one can probably construct slow blow-up 
solutions also for n < 4. 

In this section we look at radial self-similar singular solutions of i|3.1[l . i.e., the 
solutions of a special form 

V(x,t) = , , 1 rrr, W ( ; ^ | , (3.2) 

where k > 0, T £ M. are parameters. The notation of this section is not necessarily 
connected to the notation in the previous sections where we analyzed the equation 
(jl . 1|) . Substituting 1)3. 2|1 into (|3.1(l we obtain the equation for self-similar profile 

id" _)_ H—l-^yj 1 — nrw' + w 2<yJr ^ — —w — , in(0,oo). (3.3) 
r a 

Natural boundary conditions are 

«;(Q)=a, io'(0) = 0, and (3.4) 
w(r) r~j r~ x l a as r — > oo. (3-5) 

The equation 1)3. 3|l has been studied by many authors and we recommend the reader 
to consult EI El M El- 0m - aim in this section is to explain what happens to 
the solutions constructed in El an d |H1 El as we approach the critical exponent 
<ri(n) = 3/(n — 10) due to Lepin. 

We use a similar approach as we applied to (|l.l|l in the previous sections. We 
shall denote w a the solution of the initial- value problem 1)3. 3JI . 13.4JI for k — 1. We 
also define R a — sup{0 < r < oo | w a {r) > , < r < oo}. Furthermore, we introduce 
the following parameters 

1\^ If 1 

- , and /Ver = — [n — 2 

a J a \ u 

As in the previous section it is convenient to use new variables (s, it) which are defined 
by the transformation 

w(r) — r^ 1 ^ a u(r) , r — e s . 

In these variables we obtain 

^ u '-^(n-2-^)u + u 2a+l = Ke 2s u' . (3.6) 
The linearization of (|3.6(l at the trivial equilibrium u — j3 n _ a is 

~ z' + 2[n-2--)z = ne 2s z' 



a J \ a 



For the sake of brevity we denote A = n — 2 — 2/a and B = 2(n — 2 — 1/er). We 
denote Ai, A2 the roots of the characteristic polynomial for the linear second-order 
differential operator on the left-hand side. We note that the condition A > is 
equivalent to a < a s (n) = 2/ (n - 2). It is known (see @]) that l|3.3|l . I|3.4 |l -I|3.5 |l has 
no bounded solutions on (0, 00) in this case. We shall therefore assume A < (and 
hence ReAj < 0) in what follows. We also note that a > er c (n) corresponds to the 
requirement that Xj be real. 
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We define U(a, n, s) as the solution of the boundary-value problem given by Ij3.6l) 
and 1|3. 4(1 - 13. 5JI rewritten in the variables (s,u). A formal calculation, similar to the 
one leading to l|2.2L)[l in the previous section, now gives 

U(e 2s 'a, K ,s)=[3 n -Re (V lS 'e AlS Af 1 + Ai + ^, \^ S )) + ■ ■ ■ > ( 3 - 7 ) 

as s' — > oo for a suitable 7 S C, which is real if Ai is real. We emphasize that we 
always assume Re Ax < 0. 

The linearization Z at w &a solves the boundary-value problem 
n — 1 

Z" + Z' -rZ' + 2Z = 

r 

Z(0) = 1, Z'(0) = 0. 

This equation has the solution Z(r) = M(—l,n/2,r 2 /2) = 1 — r 2 /n. 
We can introduce again the index i(a) G Z 

i(a) = #{re (0, 00) I w a (r) = a a } , 

and jumps of i(a) can be used to locate the desired solutions numerically. Numerical 
computations show that i(a) jumps between 1 and 3. Another natural index in this 
case (used in PIS1I51IT2]) is 

j(a) = # (r € (0, ex.) | w a (r) = = # {s € R | U(a, 1, s) = f3 n , a } 



nl/(T 

One has j(a) = 2 for a close to a CT . The approximation (|2.20() together with some 
heuristic arguments similar to those used for justification of l|2.27(l suggest that 

when o~ s (n) < a < a c (n), 
lim j(a) = i v(n,a) when a > <r c (n) and v(n,o~) is even, 




when a > cr c (n) and v(n, a) is odd, 



where v(n, a) is defined as 



v(n } a) = min ^/cGZl ^-<fc 

This asymptotic behaviour of j(a) leads to infinitely many solutions for a s (n) < a < 
a c (n) (see [SJED- When R a < 00, j(a) must be even, hence typically we see j(a) 
jump by 2. The solutions associated with these jumps have even index j(a). 

However, the work of Lepin ([§]) shows that there are also solutions with odd 
index j(a). Our numerical calculations suggest that such solutions can be detected by 
jumps in the index i(a). The index i(a) takes on values 1 or 3. The jump 1 — > 3 of the 
index i(a) when a increases corresponds to the solutions with an even index j(cn). On 
the other hand, the jump 3 — > 1 of i(a) indicates a solution with an odd index j(a). 
These solutions cannot be numerically detected from the behaviour of j(a) only. 

Lepin showed that solutions with the index j(a) = 2 exist in the region a > a s (n), 
v{n, a) > 3. He also studied solutions with higher indices, for which he established 
existence in smaller regions of the parameter space (n,a). He conjectured that 
there are no solutions to l|3.3|l . (|3.4|l - l|3.5|l when v(n,a) < 2, which corresponds to 
n > 10 + 3/cr. 

Our numerical computations strongly support this conjecture, at least in the sense 
that no solutions from the region n < 10+3/cr can be continued outside the region. We 
looked at the solution with j(a) = 2 which appears to be the solution of l|3.3|l . I|3.4|) - 
l|3.5(l with the smallest possible a. Numerically, it also is the most robust solution. 
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Assume (at, n t , at), t G [0, 1) is a smooth path in the parameter space such that w at 
solves ()3.3|) . (I3.4() - (I3.5|1 for n = n t , a = at and K = 1. We consider the case j(at) = 2. 
and follow the path of solutions with n t — 10 — 3/a t < such that n f — 10 — 3/a t — > 0, 
nt — ► Hi, cr t — > (ji as i — > 1. Numerically, we observe that at — > -f-oo. Moreover, for 
n — 10 — 3/er>0 one expects from l|3.7|) that i? Q < +oo for large a, and we have not 
detected any solutions for numerically accessible a's. Therefore our conclusion is that 
the numerical computations provide strong evidence in favour of Lepin's conjecture. 
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